Nuclear Structure Calculations and Modern Nucleon-Nucleon Potentials 
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We study ground-state properties of the doubly magic nuclei 4 He, ie O, and 40 Ca employing 
the Goldstone expansion and using as input four different high-quality nucleon-nucleon (NN) po- 
tentials. The short-range repulsion of these potentials is renormalized by constructing a smooth 
low-momentum potential V[ ow _k- This is used directly in a Hartree-Fock approach and corrections 
up to third order in the Goldstone expansion are evaluated. Comparison of the results shows that 
they are only slightly dependent on the choice of the TVTV potential. 

PACS numbers: 21.30.Fe, 21.60.Jz, 21.10.Dr 



I. INTRODUCTION 

In recent years, the study of the properties of nuclear 
systems starting from a free nucleon-nucleon (TV TV) po- 
tential Vjvjv nas become a subject of special interest. This 
has been stimulated by the substantial progress made 
during the last decade in the development of TVTV poten- 
tials that reproduce with high precision the NN scat- 
tering data and deuteron properties 0, 0, El EI How- 
ever, the fact that these potentials predict almost iden- 
tical phase shifts does not imply, owing to their different 
off-shell behavior, that they should give the same results 
when employed in nuclear many-body calculations. It is 
therefore of great interest to investigate how much nu- 
clear structure results depend on the NN potential one 
starts with, namely to try to assess the relevance of the 
off-shell effects in microscopic nuclear structure calcula- 
tions. 

The differences between various NN potentials in de- 
scribing properties of nuclear matter have been investi- 
gated by several authors, the main aim being to try to 
assess the role of the various components of the nuclear 
force. In this context, we may mention the studies of 
Refs. 0, la 0, H i where attention has been focused on 
modern phase-shift equivalent NN potentials. 

As regards the study of finite nuclei, while a rather 
large number of realistic nuclear structure calculations 
have been carried out in the pas t few years, only a few 
attempts have been made 0, TToL Hlj to study to which 
extent these calculations depend on the NN potential 
used as input. Actually, no detailed comparison of the 
results produced by different phase-shift equivalent po- 
tentials in the description of nuclear structure properties 
has yet been done. It may be worth mentioning, however, 
the work of Ref . ^2] > where different high-precision TV TV 
potentials have been considered to investigate the effects 
of nonlocalities on the triton binding energy. 

On the above grounds, we have found it interesting 
and timely to perform nuclear structure calculations us- 
ing different phase-shift equivalent TVTV potentials and 



make a detailed comparison between the corresponding 
results. In a recent paper [T^], we performed realistic 
calculations of the ground-state properties of some dou- 
bly magic nuclei within the framework of the Goldstone 
expansion approach, and showed that the rate of con- 
vergence is very satisfactory. Motivated by the results 
obtained in that work, in the present paper we make use 
of the Goldstone expansion to calculate the binding en- 
ergies and rms charge radii of 4 He, 16 0, and 40 Ca for 
different phase-shift equivalent TVTV potentials. We feel 
that this is a good "laboratory" for a comparative study 
of the effects of TVTV potentials in finite nuclei. We con- 
sider the four high-quality TVTV potentials Nijmegen II 
0, Argonne V18 @, CD-Bonn ||, and N 3 LO ||. 

As is well known, to perform nuclear structure cal- 
culations with realistic TVTV potentials one has to deal 
with the strong repulsive behavior of such potentials in 
the high-momentum regime. Recently, a new method to 
renormalize the bare TVTV interaction has been proposed 
[tH Il5| , which is proving to be an advantageous alterna- 
tive to the use of the Brueckner G matrix |l4 . ll6ul7l , ll8j . 
It consists in deriving from Vmn a low- momentum poten- 
tial Viow-k defined within a cutoff momentum A. This is 
a smooth potential which preserves exactly the on-shell 
properties of the original Vnn ancl is suitable for being 
used directly in nuclear structure calculations. 

As in our earlier work p|, we construct the Vi ow _k for 
each of the four above mentioned TVTV potentials. The 
various Vi ow _k's are then used directly in Hartree-Fock 
calculations. Once the self-consistent basis is obtained, 
we calculate the Goldstone expansion including diagrams 
up to third order in Viow-k- 

It is worth emphasizing that one of the main advan- 
tages of the Viow-k renormalization method, with respect 
to the G-matrix one, is to preserve the phase-shift equiv- 
alence of the original TVTV potentials. Thus, it is par- 
ticularly interesting to compare the results obtained in 
nuclear structure calculations employing on-shell equiva- 
lent Viow-k 's. 

The paper is organized as follows. In Sec. II we give a 
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brief description of the main features of the four phase- 
shift equivalent NN potentials considered in our study. 
In Sec. Ill we give an outline of the derivation of Viow-k 
and some details of our calculations. In Sec. IV we 
present and discuss our results. Some concluding remarks 
are given in Sec. V. 



II. REALISTIC NUCLEON-NUCLEON 
POTENTIALS 

As reported in the Introduction, we employ the Ni- 
jmegen II pj, Argonne V18 0, CD-Bonn and chiral 
N 3 LO NN potentials, which have all been fitted to 
the Nijmegen phase-shift analysis as well as the proton- 
proton and neutron-proton data below 350 MeV [l9j . 
It is well known that these potentials, even if they re- 
produce the NN data with almost the same accuracy, 
may have a rather different mathematical structure. The 
Nijmegen II and the Argonne V18 potentials are non- 
relativistic and defined in terms of local functions, which 
are multiplied by a set of spin, isospin and angular mo- 
mentum operators. The CD-Bonn potential, based on 
relativistic meson field theory, is represented in terms 
of the covariant Feynmann amplitudes for one-boson ex- 
change, which are nonlocal [2C|. The N 3 LO potential 
is based upon a chiral effective Lagrangian. The model 
includes one- and two-pion exchange contributions and 
so-called contact terms up to chiral order four, some of 
which are non-local. 

All the above NN interactions reproduce equally well the 
same phase-shifts, so the corresponding on-shell matrix 
elements of the reaction matrix T are the same as well. 
However, this does not imply that the interactions are 
identical. The T-matrix is obtained from the Lippmann- 
Schwinger equation 



T(k',k,k 2 ) = V NN {k',k) + 



-V 



q 2 dqV NN {k',q)- 



T T(q,k,k 2 ) , (1) 



where k, k', and q stand for the relative momenta. No- 
tice that the T-matrix is the sum of two terms, the Born 
term and an integral term. Notwithstanding the sum is 
the same, the individual terms may still be quite dif- 
ferent. For example, in Fig. 1 we show the value of 
the 3 Si T-matrix elements for K\ a \> = 150 MeV and 
k = k' = ko = 1.34 fin - . They are indicated by the 
full circle and are practically identical for the four NN 
potentials we have employed. In Fig. 1 we display also 
the matrix elements Vivjv(&o,fc) as a function of k for 
the four potentials. The asterisks stand for the diagonal 
matrix elements Vjvj\r(M)! &o)j which represent the Born 
approximation to T. From the inspection of this figure, 
it is clear that even if different Vnn's reproduce the same 



T-matrix element, the latter is obtained by summing two 
terms that are significantly different for each potential. 
More precisely, it is evident that Nijmegen II potential 
is a quite "hard" potential, its diagonal matrix element 
being very repulsive. So, in order to reproduce correctly 
the on-shell T-matrix element, it needs a large attractive 
contribution from the integral term of Eq. (1), the latter 
being related to the tensor component of the NN force 
and to its off-shell behavior (see Ref. |2l| for a closer ex- 
amination). On the other hand, the N 3 LO interaction is 
a rather "soft" potential, implying a smaller contribution 
from the integral term. 
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FIG. 1: Matrix elements VjvAr(fco,fc) for the 3 Si partial wave 
for CD-Bonn (green line), Nijmegen II (red line), Argonne 
V18 (black line), and N 3 LO (blue line) potentials. The di- 
agonal matrix elements ko = 1.34 frn -1 are marked by an 
asterisk. The corresponding T-matrix element is marked by 
a full circle (see text for more details). 



Similar features may be observed also in the other par- 
tial waves, however, the differences among the potentials 
decrease for larger values of the orbital angular momen- 
tum. 



III. METHOD OF CALCULATION 

A traditional approach to the renormalization of the 
strong repulsive behavior of realistic NN potentials, 
when dealing with doubly closed-shell nuclei, is the 
Brueckner-Goldstone (BG) theory (see for instance Refs. 
[22I Hj|), where the Goldstone perturbative expansion 
is re-ordered summing to all orders only the ladder dia- 
grams. Consequently, the bare interaction (Vnn) vertices 
are replaced by the reaction matrix (G). This framework 
leads to the well known Brueckner-Hartree-Fock (BHF) 
theory, when the self-consistent definition is adopted for 
the single-particle (SP) auxiliary potential and only the 
first-order contribution in the BG expansion is taken into 
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account. So, the BHF approximation gives a mean field 
description of the ground state of nuclei in terms of the 
G matrix, the latter taking into account the correlations 
between pairs of nucleons. However this procedure is not 
without difficulties, because of the energy dependence of 
G. 

As already mentioned in the Introduction, we renor- 
malize the short-range repulsion of the bare NN poten- 
tial by integrating out its high momentum components 
|14lll5l |. The resulting low-momentum potential, Vjow-k, 
is a smooth potential, whose vertices can be used directly 
to sum up the Goldstone expansion diagrams. 

According to the general definition of a renormaliza- 
tion group transformation, Vi ow -k must be such that the 
low-energy observables calculated in the full theory are 
preserved exactly by the effective theory. 

For the nucleon-nucleon problem in vacuum, we re- 
quire that the deuteron binding energy, low-energy phase 
shifts, and low-momentum half-on-shell T matrix calcu- 
lated from Vnn must be reproduced by Viow-k- The 
effective low-momentum T matrix is defined by 

Tiow-k(p',P,.p 2 ) = Vi ow -- k {p',p) + 



+V / q 2 dqVi ow - k (p',q)- 



T Ti ow _ k (q,p,p 2 ) . (2) 



Note that for Xi ow -k the intermediate states are inte- 
grated up to A. 

It is required that, for p and p' both belonging to 
PfEitf < A), T(p',p,p 2 ) = T low _ k (p\p,p 2 ). In Refs. 
[pa, Il5| it has been shown that the above requirements 
are satisfied when Vj ow -k is given by the folded-diagram 
series 



ow-k 



= Q~Q 



■JQ + Q'J 



Q / Q 



TABLE I: Deuteron binding energy (MeV) and np So phase 
shifts (deg) as predicted by full CD-Bonn and its Vi ow _k (A = 
2.0 fin -1 ) 
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2.224 


2.224 


2.224 
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62.1 


62.1 


62.1 


10 


60.0 


60.0 


60.0 


25 


50.9 


50.9 


50.9 


50 


40.5 


40.5 


40.5 


100 


26.4 


26.4 


26.8 


150 


16.3 


16.3 


16.9 


200 


8.3 


8.3 


8.9 


250 


1.6 


1.6 


2.0 


300 


-4.3 


-4.3 


-4.5 



with the full CD-Bonn potential and its Vi w-k (with a 
cut-off momentum A = 2.0 fin -1 ). 

An important question in this approach is what value 
one should use for the cutoff momentum. A discussion of 
this point as well as a criterion for the choice of A can be 
found in Ref. [l5|- According to this criterion, we have 
used here A = 2.1 fm~ . 

After having renormalized the various NN potentials, 
we use the corresponding Vi ow _k's directly in a HF cal- 
culation. The HF equations are then solved for 4 He, 16 O 
and 40 Ca making use of a harmonic-oscillator (HO) ba- 
sis. The details of the HF procedure are reported in Ref. 
|13|. As a major improvement, in this work we remove 
the spurious center-of-mass kinetic energy writing the ki- 
netic energy operator T as 



T 



2Am 



(4) 



Similarly, we define the mean square radius operator 



-Q'jQ I Qj Q + ... , (3) 

where Q is an irreducible vertex function, in the sense 
that its intermediate states must be outside the model 
space P. The integral sign represents a generalized fold- 
ing operation [24| , and Q' is obtained from Q by removing 
terms of first order in the interaction. 

The above Vi ow _k can be calculated by means of itera- 
tive techniques. We have used here an iteration method 
proposed in Ref. p5|. which is particularly suitable for 
non-degenerate model spaces. This method, which we 
refer to as Andreozzi-Lee-Suzuki (ALS) method, is an 
iterative method of the Lee-Suzuki type p6| . 

To exemplify how Vi w-k preserves low-energy observ- 
ables, we report in Table I the deuteron binding energy, 
and the neutron-proton So phase shifts calculated both 



A complete review about center-of-mass correction in 
self-consistent theories may be found in Ref. j2^| 

In our calculations the HF SP states are expanded 
in a finite series of N — 5 harmonic-oscillator wave- 
functions for 16 and 40 Ca, and N = 6 for 4 Hc. This 
truncation is sufficient to ensure that the HF results do 
not significantly depend on the variation of the oscilla- 
tor constant fkv, as we showed in Ref. ^3|. The values 
of hid adopted here have been derived from the expres- 
sion Two = 45A~ 1 / 3 - 2hA~ 2 / 3 28], which reproduces 
the rms radii in an independent-particle approximation 
with harmonic-oscillator wave functions. This expression 
gives hw = 18, 14 and 11 MeV for 4 He, 16 and 40 Ca, 
respectively. 
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We use the HF basis to sum both the Goldstone expan- 
sion and the diagrams for the mean square charge radius 
(r 2 ), including contributions up to third order in Viow-k- 
Fig . 2 shows first-, second-, and third-order diagrams 
|29| of the Goldstone expansion. 




(c) (d) (e) 



FIG. 2: First-, second-, and third-order diagrams in the Gold- 
stone expansion. 



IV. RESULTS 

In Table II we show for 4 He, 16 and 40 Ca the calcu- 
lated binding energy per nucleon and the rms charge ra- 
dius obtained using different phase-shift equivalent NN 
potentials, and compare them with the experimental data 



TABLE II: Comparison of the calculated binding energies per 
nucleon (MeV/nucleon) and rms radii (fm) for different Vnn 
with the experimental data for 4 He, le O and 40 Ca. We take 
into account the finite dimensions of the proton using an es- 
timate of its rms charge radius \fr^ = 0.8 fm |3ll ]. 



Nucleus Nijmegen II AV18 CD-Bonn N3LO Expt. 



BE /A 

(r 2 ) 1 ' 2 


6.88 


6.85 


6.95 


6.61 


7.07 


1.68 


1.69 


1.63 


1.75 


1.67 


16 












BE /A 

(r y/2 


8.26 


8.26 


8.30 


8.11 


7.98 


2.58 


2.59 


2.49 


2.66 


2.73 


40 Ca 












BE /A 

(t 2 ) 1 ' 2 


9.66 


9.53 


9.93 


9.50 


8.55 


3.20 


3.22 


3.10 


3.29 


3.485 



A detailed analysis about the convergence properties of 
the perturbative series can be found in Ref. \\?\ , where it 
is shown that the convergence is fairly rapid and higher- 
order contributions are negligible. 

From Table II, we see that the calculated quantities 
are scarcely sensitive to the choice of the NN potential. 

As a matter of fact, the binding energies per nucleon 
and the rms charge radii calculated using the various po- 
tentials differ at most by 0.43 MeV and 0.19 fm, respec- 
tively. This insensitivity may be traced back to the fact 
that when renormalizing the short-range repulsion of the 



TABLE III: Calculated P D 's with different NN potentials 
and with the corresponding Vj ow _k's 





Nijmegen II 


AV18 


CD-Bonn 


N3LO 


Full potential 


5.63 


5.76 


4.85 


4.51 


Fiow-k 


4.32 


4.37 


4.04 


4.32 



various potentials, the differences existing between their 
off-shell properties are attenuated. It is well known (see 
for instance [2l]]) that the off-shell behavior of a poten- 
tial, and in particular its off-shell tensor force strength, 
is related to the Z?-state probability of the deuteron Pjy; 
this is why when comparing NN potentials, off-shell dif- 
ferences are seen in Pjj differences. For this reason, we 
report in Table III the predicted Pd's for each of the 
potentials under consideration, and compare them with 
those calculated with the corresponding Vjow-k's. We 
see that while the Pd's given by the full potentials are 
substantially different, ranging from 4.5 to 5.8 %, they 
become quite similar after renormalization. This is an 
indication that the "on-shell equivalent" potentials we 
have used are made almost "off-shell equivalent" by the 
renormalization procedure. This aspect is more evident 
when comparing our calculated 4 He binding energies with 
the results of exact calculations, based on the Faddeev- 
Yakubovsky procedure [^3, 0] (see Table IV) . As a mat- 
ter of fact, the difference between the exact results, that 
is at most 2 MeV, is reduced to 1.3 MeV in our calcu- 
lations. It has to be observed that larger differences be- 
tween our results and the exact ones are obtained when 
employing high-P^i potentials, such as Nijmegen II and 
Argonne V18. Smaller differences occur for the N 3 LO (1 
MeV) and CD-Bonn (1.5 MeV) potentials, whose tensor 
force strengths are smaller than those of the two other 
potentials. This reflects the fact that for these two poten- 
tials the renormalization procedure modifies the original 
Pd to a limited extent, while a stronger change occurs 
for the Nijmegen II and Argonne V18 potentials. Table 
IV also shows that in all cases we get more binding than 
the exact calculations, as a consequence of the renormal- 
ization of the repulsive components of the potentials. 



TABLE IV: Comparison of the 4 He calculated binding ener- 
gies (MeV) for different Fnn obtained using Goldstone expan- 
sion (I) and with Faddeev-Yakubovsky procedure (II) . 



BE 


Nijmegen II 


AV18 


CD-Bonn 


N3LO 


(I) 


27.523 


27.409 


27.799 


26.440 


(II) 


24.560 


24.280 


26.260 


25.410 



V. SUMMARY AND CONCLUSIONS 

The aim of this work has been to compare the results of 
microscopic nuclear structure calculations, starting from 
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four different phase-shift equivalent NN potentials, Ni- 
jmegen II, Argonne V18, CD-Bonn, and N 3 LO. To this 
end, we have calculated ground-state properties of the 
doubly closed nuclei 4 He, 16 O, and 40 Ca by way of the 
Goldstone expansion. This has been done within the 
framework of the so-called Mow-k approach [hH ITfij to 
the renormalization of the short-range repulsion of the 
NN potentials, wherein a low-momentum potential is 
derived, which preserves the low-energy physics of the 
original potential. 

The analysis of the results obtained shows that the 
calculated properties are only weakly dependent on the 
NN potential used as input. This result may be traced 
back to the renormalization procedure of the short-range 
repulsion. As a matter of fact, we have shown that the 
renormalized potentials are characterized by a reduced 
off-shell tensor force strength, as compared with that of 
the original potential. Moreover, the Pd's of the differ- 
ent potentials become quite similar, this quantity being 
related to the balance between the central and tensor 



components of the nuclear force. 

It is worthwhile to point out that when dealing with 
the N 3 LO chiral potential, the renormalization proce- 
dure, which by construction preserves exactly the on- 
shell properties up to the cutoff momentum, scarcely 
modifies the off-shell behavior. This feature, which is 
related to the fact that chiral perturbation theory is a 
low-momentum expansion, may make this kind of NN 
potentials particularly tailored for microscopic nuclear 
structure calculations. 
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